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AND 
K. Y. LO 


Norwegian Geotechnical Institute, Oslo 


SUMMARY 


A theory of one-dimensional consolidation taking account of secondary time 
effects is presented. General expressions for excess pore water pressure and settle- 
ment for any type of monotonic time dependent loading have been obtained and 
these are specialized to the case of step loading. 

The parameters involved in the theory are the permeability, the viscosity, and 
the primary and secondary compressibility of the soil structure. Methods for deter- 
mining these parameters from laboratory consolidation tests are discussed. 

Some laboratory test data on three clays: London clay, Grangemouth clay, 
and Sodium Bentonite are presented which seem to justify the assumptions made 


in the theory. Finally, the application of the theory to settlement analysis is con- 
sidered. 


he 


NOTATIONS 


(L= Length, T= Time, W = Weight, 0 


compressibility of spring ‘a’ 


(primary compressibility) [L*/W] 
) id 
viscosity of dashpot ‘A’ (WT/L?*) 
vertical effective stress 
strain in ‘b’ and ‘A’ [0] 
time clapsed after application of (T] 


pressure 


applied pressure, a function of time [W/L?*] 


excess pore water pressure [W/L?*]} 
total pressure [W/L] 
permeability of soil {L/T] 
density of water [W/L] 
i ral transfor ra- 

integral transform para 
Laplace transform of o’(t) (WT/L?] 
Laplace transform of q(t) [WT/L?} 


Xn 
U 
e(t) 


dimensionless.) 


coefficient of consolidation 


A((1/a) + (1/b)] {1/T] 
Alb {1/T] 
[p(p + a)]/[O(p + [T/L] 
[(O/h*) (a®/4)] n® where n is 
any positive odd integer [1/T] 
[a/(a + b)) K, ({1/T) 
(a+ + K,)*-4pK,) (1/7) 
compressibility factor = 1 + (b/a) [0] 
influence factor = Ah? /bO [0] 
primary time factor = Ot/h® [0] 
(4h?/ 0) [0] 
settlement at any time ¢ {L] 
average degree of consolidation [0] 
compressive strain at any time / [0] 


e(%) compressive strain at infinite time [0] 
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1, INTRODUCTION 


In the classical theory of one-dimensional consolidation developed by 
Trerzacut (1923) the relation between the void ratio and the vertical effective stress 
acting in the soil skeleton is assumed to be intrinsically time independent. This 
theory leads to results substantially in accord with laboratory test data for a wide 
range of clays, but for certain soils important departures have been observed (see, 
for example, Taylor 1942). In particular, it has been found that the rate of settle- 
ment decreases toward the end of the consolidation process less rapidly than would 
be expected from the theory, and this additional settlement has been termed 
“secondary 

Numerous attempts have been made to provide a basis for predicting, from 
the results of laboratory consolidation tests, the magnitude of secondary settlement 
likely to occur in loaded natural clay strata, and these approaches fall into two 
groups: the empirical and the analytical. Both laboratory tests and field settlement 
records (BuisMAN 1936, Taytor 1942, and Scuaap 1948, Kopresan 1948, 
Gruze 1957, Zexvaert 1958) suggest that secondary settlement becomes important 
only when the excess pore water pressures have practically dissipated, and that 
it increases proportionally with the logarithm of time. These observations form the 
basis of some early empirical methods (Buisman 1936, Koppejan 1948) and 
although they are valuable for certain purposes and have the advantage of being 
readily applied in practice—a virtue which the more elaborate analytical theories 
do not possess—they all imply that secondary settlement is non-terminating. More- 
over, the time at which secondary consolidation commences is not clearly defined. 
Later attempts were made, notably by Zeevaert (1958) and SuKtyje (1957), to place 
these empirical methods on a more rational foundation; these theories may be 
termed semi-empirical. 

The first rational theories were formulated by Taytor and Mercuanrt (1940), 
Taylor (1942), and later by Goldstein (1952). 

More recently important contributions have been made by Tan 1957 (a, b, c, d). 
McNass (1960) has derived the equation governing one-dimensional consolidation 
in a very general form which covers the case of finite compression and permits the 
natural incorporation into the analysis of such effects as varying compressibility, 
permeability and creep. In these theories the basic ideas underlying Terzaghi’s 
theory are retained, but in addition it is supposed that the soil skeleton exhibits 
creep under constant effective stress. While the consequences of assuming the soil 
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observed settlement characteristics of clay specimens exhibiting 
secondary effects, the difficulty of determining the relevant soil 
parameters from the test data has limited the practical application of 
these ideas. Even more involved are the three-dimensional theories 
advanced by Biot (1956), MANpeEL (1957), and TaN (1957 e¢, f). 
In this paper we shall be concerned only with one-dimensional 
consolidation and we shall suppose that the soil skeleton can be 
o'() replaced by a certain three clement rheological model which was 
originally suggested by Tan (1957 a); it is equivalent to a special 
Fig. 1. Model case of Biot’s (1956) general theory. 
of sail diddoen. The equations governing one-dimensional consolidation in such 
an idealized soil will be formulated, and expressions derived for the 
settlement and the excess pore water pressure distribution at any time. Methods for 
determining the relevant soil parameters are suggested and typical experimental 
results are presented, but we defer a full discussion of experimental results to a later 
paper. 


{ skeleton to behave in this way are in broad agreement with the 


2. SCHEMATIC REPRESENTATION OF THE CONSOLIDATION PROCESS 


The proposed model of the soil skeleton undergoing simple one-dimensional 
compression consists of a Hookean spring connected in series with a Kelvin or Voigt 
element as shown in Fig. (1). We first examine how the physics of consolidation 
can be correlated with the behaviour of this simple model. When a time dependent 
stress o’(¢) acts on an element of the soil at a depth z, where the response of the soil 
skeleton is schematically represented by the model, the spring ‘a’ compresses 
instantaneously, but the settlement of the compound element ‘b’—‘A’ is retarded 
because of the presence of the dashpot ‘A’. Since the transference of stress from 
pore water to soil skeleton is delayed due to the low permeability of the clay, 
the effective stress o’(¢) increases gradually from zero to the full value of the 
applied stress. Hence, the compression of spring ‘a’ is also gradual and is fully 
accomplished only when the applied stress has become fully effective. 

The period during which the settlement and rate of settlement is dominated 
by compression of spring ‘a’ is often called the hydrodynamic period and the 
consolidation that occurs the primary consolidation. 

Under the gradually increasing effective stress the Kelvin body commences 
to compress. Initially the full load is taken by the Newtonian dashpot ‘X’ and is 
progressively transferred to the spring ‘b’ as compression proceeds. This pheno- 
menon of transference corresponds to the process of secondary consolidation which 
occurs under sustained effective stress. After a very appreciable time has elapsed 
the full effective stress will be taken through springs ‘a’ and ‘b’ the dashpot ‘A’ 
sustaining no load. 
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38. FORMULATION OF PROBLEM 


At any time ¢ after the application of a load, the effective stress acting on the 
model will be o’(t). The compressive strains in the elements ‘b’ and ‘A’ must be 
equal, The strain f in ‘b’ is proportional to its load while in ‘A’ the load governs 
the rate of viscous deformation. Since the sum of the stress in elements ‘b’ and‘ dA’ 
must be equal to the effective stress o’(t) it follows that 


1d 


In the above equation both f and o’(¢) are functions of the depth z. The solution 
of this equation, satisfying the initial condition f(0) = 0 is 


(1) 
fm Pas (2) 
0 
The effective stress-strain relationship can therefore be expressed in the form 
e=a0'() + (3) 
0 


If we consider the continuity of flow of water through an element of soil (Fig. 2) 
then it can be shown that for small strains 


The pore water pressure p’ and the effective stress 0’ are related to the applied 
stress 0 by the well-known equation 
=g- p’ 
Moreover the excess pore water pressure may be defined by 
p = - 20) 
where the sign is chosen according as z is measured upward or downward and 
z, denotes the position of the water level. 


Since one-dimensional consolidation is strictly compatible only with an applied 
pressure change o varying linearly with depth it follows that 


k de (6) 


ot 


q (t) 
z 


Fig. 2. 


of a finite layer of clay. impermeable base oz 
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Combining equations (3) and (6) we obtain 


k d%a’ da’ 
0 


which is the one-dimensional equation of consolidation expressed in terms of the 
effective stress. 


This equation, together with the boundary conditions 


0O<tS o (8) 
= q(t), 0 


completely defines the problem under consideration. 


4. SOLUTION FOR MONOTONIC LOADING 


We consider a general monotonic loading programme in which the applied 
surface pressure is an arbitrary function of time q(t). The initial pressure (0) is 
applied at the instant ¢= 0. We denote the Laplace Transform of 0'(z,t) by o'(z,p), 
that is 

t 


a’ (z,p) = Jo (z,t)e (9) 


0 


Applying a Laplace Transform to both sides of equation (7), and using the 


Convolution Theorem for the last term on the right hand side, the ordinary 
differential equation 


d?a’ 
ie =0 (10) 
is obtained 
_ p(p + 
where O(p+8)’ 
(- 5) 
p= 
aYw 
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The general solution of (10) is 
= A sinh rz + B cosh rz (11) 


where A and B are arbitrary functions of p, which must be determined from the 


transformed boundary conditions (8) which are 


da’ 
= 4p), 
It follows that 
A = ~q(p) tanh rh, 
B= q(p), 
h z) 
at rh (13) 


and hence 
which is the formal solution for the vertical effective stress distribution through the 


thickness of the layer. 


5. EXPRESSION FOR SETTLEMENT AND PORE WATER PRESSURI 
In practice, the quantities which are of interest to the engineer are the 


settlement S(t) of the surface of the layer at any time ¢, and the pore water 


pressure p’(z,0). 
he settlement S can be obtained by transforming and integrating (3) 
S(p) = J e(zp) de (14) 


which gives 
= alo p+ 
S(p) = 5) q(p) tanh rh (15) 
The settlement may then be represented by the contour integral 
ay® 


and the pore water pressure is given by 


cosh r (h—-z) 
(x,t) = q(t) | coshth (17) 
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Expressions (16) and (17) are valid for all types of monotonic loading, and 
when q(t) has been specified, explicit solutions for the settlement and the pore 
water pressure may be found. We shall now consider the simplest type of loading: 
a step load such as shown in Fig. (3). 


The loading function is specified by 
q(t) = 0 t<0. 
q(t) = qo t>0. 


with the corresponding transform 


q(b) = (18) 
From (15) and (18) it follows that 
19 
S(p) = (2) tanh rh - (19) 
Using the known expansion 
= 
tanh x =8x > (20) 
n= odd 
equation (19) then becomes 
© 
n= odd 


The image of the expression under the summation sign is 


Fe 
where + Ki) + KC + K,)*—46K;] (22) 
XQ 2 
x? @ 
= (23) 


The complete inversion of (21) from the p-plane to the ¢-plane can now be 
effected and the final expression for the settlement is found to be: 


X4— XQ XQ 
n= odd 


a 
where K = 


and the pore water pressure is found to be given by 


8 
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Fig. 3. Loading programme. 
q(t) 

4, 

t 
1 f4atb 1 “a+b 
n — Xp 
n= odd 
a x 

(- +b — set , (25) 
sin 

Xy — 2h 


For soil which is infinitely viscous (A > 0) or when there is no secondary com- 
pression (b > 0), the above equations reduce to 


D 


S(t) —aqh 4 (24 a) 
n= odd 
and 
D 
4 Ot 95 
p te 4 sin wh ( ) 
n= odd 


These may be identified as the classical solution given by Terzaghi; the results 
might have been anticipated from purely physical reasoning. 


6. APPLICATION OF THEORY 


The principal application of consolidation theory is the prediction of the 
amount and time rate of settlement of engineering structures. For this purpose, 
it is necessary to 


(a) determine by suitable laboratory experiments the magnitudes of all the para- 
meters involved in the theory. 

(b) choose the appropriate theoretical time-settlement curve appropriate to the 
field case under consideration. 


The theory involves four parameters: the permeability k, the viscosity of the 
soil structure 1/A, the primary compressibility a and the secondary compressibility 
b. To determine these parameters, which describe the consolidation characteristics 
of a soil, we shall need expressions for S(t) appropriate for large and small values 
of time. These may be obtained by suitable recasting and inversion of equation (15). 
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The solution for small values of time can “ shown to be 


S(t) = [14+ 5% —t (26) 


and for large values of time, the following expression has been formed 


Si) ¢ | (27) 
We observe that in the limit / > 00, 


(a+ (28) 


as it should. 


Experimental evidence has shown that for remoulded London clay (L.L.= 75"/o, 
P.L. = 30%), Grangemouth clay (L.L. 62 °/o, P.L. = 26 °/o), and Sodium Bentonite 
(L.L. = 530 %o, P.L. = 120 °/o), the value of A/a is very small compared with unity 
(approximately 5  10-° 1/min). For other soils, the same order of magnitude may 
be expected. Hence, the coefficient of consolidation ©[= k/(ay,)| can most easily 
be determined by plotting the settlement against the square root of time; this method, 
of course, is used conventionally. 

From equation (27), the parameters a, b, A, may be found in the following way: 
We write (27) in the form 


=(a+b)—be (29) 
or alternatively from (28) we have 
e(co)—e(t) _ (30) 


Taking the logarithm of both sides 


logiob — 0,434 ~ (31) 


logio 

Hence a plot of log [e(#%) — e(t)]/o’ against t will yield the intercept b and the 

slope A/b. The value of a may be determined either from the final settlement or 
from the relation 


A 
a’ 

In Fig. (4) is shown a typical plot which may be expected from experimental 
data. The full line is the curve that would actually be obtained from test data 
while the dotted line is that which would have been found if there were no hydro- 
dynamic time lag. Indeed, equation (27) can be shown to be the solution for the 
case when the excess pore water pressure has virtually dissipated. This equation 
therefore holds in the range of secondary consolidation when the applied stress 
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Fig. 4. Slope of BD = 0,434 A/b 
Determination of Soil c Intercept OB = log b 
Parameters ~ Method 1, Fort *t, , u O, 
a 1 b a 
ofte? 
-Cutve trom test data 
8 Also: a b a’ 
OTT A 
0 


0 


Time t (min) 


may be considered to be fully effective. This concept is in accord with the experi- 
mental evidence obtained by Taylor (1942), Koppejan (1948) and more recently 
by Geuze (1957). These authors pointed out that secondary compression takes place 
at a time when the excess pore water pressure is practically negligible.. The 
physical significance of Fig. (4) is therefore apparent. The emergence of the curve 
into a straight line at point A(t = ¢,) gives the approximate time at which the excess 
pore water pressure may be considered to have practically dissipated. Figures (5) 
to (7) show some typical experimental data on three clays: London clay, Grange- 
mouth clay and Sodium Bentonite; these plots give not only the values of the soil 
parameters but provide an indirect verification of the theory. 

In figures (5) and (6), some curves were plotted using an alternative method 
which has the same basis as that discussed. In this method, log ,(4» — 4;) is plotted 
against ¢;Here 8, and 4, are the settlements at times /, and ¢; respectively. Using this 
method, it is not necessary to obtain the final settlement. Values of the soil para- 
meters obtained by both methods are identical within limits of experimental 
accuracy. 

It is, however, too lengthy to proceed to a detailed discussion of experimental 
results here. To use the results of our investigation to full advantage it is necessary 
to derive an expression for settlement in dimensionless form. To this end we 
introduce the following dimensionless parameters: 


M 
ah? 
= 33 
(38) 
Ot 
T, 


We shall call M the compressibility factor, N the influence factor, and T, 
the primary time factor which is identical with Terzaghi’s time factor. Equation 
(24) may then be written 
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Increment] | |ib/in®/min| i/min | Method | Temp*c 


ndon Cli rown 
Lb. * 74,7 #70,0 
PL.+29,7 @, #1,920 
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a b 


0,8 26x10? (1) 20 
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Fig. 5. Determination of Soil Parameters of London Clay. 
Grangemouth Clay 
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t(minx10*) (Method 1) 
or: t (min x10?) (Method 2) 


Fig. 6. Determination of Soil Parameters of Grangemouth Clay. 
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S(t) 
nm odd 


Xe Xa 2 


(35) 


The expression (34) gives, in effect, the degree of consolidation, so that it may 
be written 


U = F(M, N, 7;) (36) 


It is evident, therefore, that the consolidation characteristics of the ideal soil 
we have considered depends not only on the time factor 7,, but also on two 
additional parameters M and N. It should be noted that N is a function not only 
of the soil constants but also of the thickness / of the consolidating stratum. 


A set of theoretical U - 7, curves is shown in Fig. (8) for M= 1.10 and for 
a range of values of N. Hence, M and N may be determined from laboratory tests, 
taking due account of the thickness and the appropriate field consolidation settle- 
ment-time curve may be chosen. 


A further point of practical significance must be mentioned. As the thickness 
of a laboratory sample is generally much less than that of a naturally occurring 
clay stratum, the N values for these two layers will be quite different, and the 


| | Bentonite { 
Test No} 
o | w b s/o | 
\ increment Symbol % | in@/ib in2/ib- i/min | Method | Temp | 
o 3 0,10 | 844 1,0x10 (1) | 
x $42 92410) 0,91x10"? 3.02105 (1) | 
— 
2 | | | 
SS 
25 
lo 
20 
aa! 
1,0} 
0,5 
13 
T(min x 109) 
Fig. 7. Determination of Soil Parameters of Sodium Bentonite. 
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Theoretical Curves of U(M,N,1) 
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b 4 
An? 
T 
— 
+4 
10? 104 10° 


Primary Time Factor 1, 


Fig. 8. Theoretical Curves of average degree of Consolidation, 
U for a soil having a secondary compression 10 "/o that of the 


primary settlement. 


shapes of the settlement-time curves for these two cases will not be simply related. 
Only when secondary consolidation is absent (A = 0 or 0) and N assumes a definite 
value (0 or ©), independent of h, will the shapes of these two curves be simply 
related by the “ratio of the squares of the thickness” rule. 


7, CONCLUSIONS 


A one-dimensional theory of consolidation taking secondary effects into con- 
sideration has been developed. The assumptions of Terzaghi’s classical theory are 
retained with the exception that the soil skeleton now responds to changes of 
effective stress like a spring in series with a Kelvin body; the effective stress-strain 
relationship involves therefore an intrinsic time dependence. The consolidation 
phenomenon thus consists of two time dependent processes: the dissipation of excess 
pore water pressure and the creep of the soil skeleton under constant effective stress. 
The mathematical discussion establishes that the average degree of consolidation 
is related not only to the primary time factor 7), but also to two dimensionless 
parameters: the compressibility factor M and the influence factor N. These para- 
meters characterize respectively the magnitude and rate of secondary consolidation. 
Some typical U - 7, curves have been plotted for a range of values of N for the 
case when the ultimate secondary settlement amounts to 10 °/o of the primary 
settlement. 

The theoretical study also leads to an expression for the distribution of excess 
pore water pressure and settlement at any time following the load application. 
These formal expressions are valid, subject to the restriction that no unloading 
occurs, for a load varying with time in an arbitrary manner. 
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Methods have been provided whereby all the parameters involved may be 
determined from experimental data without recourse to “fitting” methods. Test 
data so far obtained from the three clays cited confirm the reasonablenesss of the 
basic assumptions made. The application of the theory to practical problems can 
easily be achieved by laboratory determination of the soil consolidation parameters 
a, b, and A, and these, together with the thickness of the layer under consideration, 
defines the values of M and N and thus enables the field settlement-time curve to 
be predicted. 
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